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Abstract 

In this survey paper we present a class of shape optimization problems where the 
T-H cost function involves the solution of a PDE of elliptic type in the unknown domain. 

lO In particular, we consider cost functions which depend on the spectrum of an elliptic 

operator and we focus on the existence of an optimal domain. The known results are 
O presented as well as a list of still open problems. Related fields as optimal partition 

problems, evolution flows, Cheeger-type problems, are also considered. 
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U 

Q 1 Introduction 

Looking for optimal shapes is a very fascinating field, perhaps due to the fact that a shape is 
^ something closer to the human spirit than a function, in which a parametrization procedure 

' is present, making it less direct and intuitive. However, from a mathematical point of view, 

handling shapes is much more difficult than handling functions, for instance because no 
I> vector space structure can be defined on a family of domains, which makes useless most of 

the functional analysis tools in Banach spaces. 
A general shape optimization problem reads as 

m 

mm{F{n) : QgA) (1.1) 

^) where ^ is a suitable family of admissible domains and F is a suitable cost function defined 

on A. Problems of this kind arise in many fields and in many applications, and the literature 
^ is very wide, from the classical cases of isoperimetric problems and the Newton problem of 

^ the best aerodynamical shape to the most recent applications to elasticity and to spectral 

^ optimization. Without any claim to be complete and exhaustive we quote the recent books 

[31 [m [261 EZl EQI Uni ESI) where the reader can find a lot of shape optimization problems 
together with all the necessary details and references. 

A family of shapes is in general more rigid than a family of functions, in the sense 
that a sequence of shapes may have a limit which is no more a shape; then quite often 



problems like (1.1) do not have optimal solutions and optimal shapes do not exist, even for 
very natural cost functions F and admissible domains A. In these situations, in order to 
perform the analysis of minimizing sequences, a relaxation procedure is needed, enlarging 
the class A to objects that are only limits of shapes and defining the cost F on these new 
objects. 
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It is not our goal to describe shape optimization problems and their relaxed formulations 
in full generality; we will limit ourselves to the cases where the cost F is of the form 

F{a) = J{un) 

being uq the solution of some elliptic PDE in 0. Two main cases fall for instance in this 
scheme: the case of integral functionals and the case of functions of the spectrum of an 
elliptic operator. 

Integral functionals. Given a right-hand side / we consider the PDE 

-Au = fmn, u£H^{n) 

which provides, for every admissible domain O C M*^, a unique solution uq that we assume 
extended by zero outside of 17. The cost F{Q) = J{u^) is obtained by taking 

J(u)= / j{x,u{x))dx 

for a suitable integrand j. 

Spectral optimization. For every admissible domain Q consider the Dirichlet Laplacian 
—A which, under mild conditions on Q, admits a compact resolvent and so a discrete 
spectrum A(0). The cost is in this case of the form 

F{n) = <^{x{n)) 

for a suitable function $. For instance, taking $(A) = we may consider the optimization 
problem for the fe-th eigenvalue of —A: 

mm{Xk{0.) : A}. 

We will see that, adding suitable geometrical restrictions to the class A of admissible 
domains, gives some extra compactness properties and prevents the relaxation of minimizing 
sequences, thus providing in many cases the existence of optimal domains. In some other 
cases, the existence of optimal domains will be obtained without any geometrical restriction, 
as a consequence of some qualitative properties of the cost function. 

The plan of the paper is as follows. 

o In Section [2] we introduce all mathematical tools that are necessary to treat the 
optimization problems in the rest of the paper. In particular capacity, capacitary 
measures, 7-convergence are recalled, together with their main properties. 

o Section [3] deals with optimization problems for integral functionals of the form above. 
We show that in general one cannot expect that an optimal domain exists and the 
optimum has to be searched in the class of capacitary measures, that is the relaxation 
of the class of domains. Some necessary conditions of optimality are also provided. 

o In Section 4 we consider the case of optimization problems for a function $(A(17)) 
where A(r2j is the spectrum of an elliptic operator, that we always take —A with 
Dirichlet conditions at the boundary. 

o Section [5] deals with the particular case of cost functions which only depends on the 
first two eigenvalues, of the form $(Ai(r2), A2(17)) . The particular form of the cost 
allows to obtain the existence of optimal domains under very mild assumptions on 
the function 
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o Section [6] is devoted to optimal partition problems, where the unknown is a partition 
Ui^i of ^ given set D in a given number N of domains, and the cost function 
F(r2i, . . . is of a fairly general type. Some interesting cases of spectral costs are 
presented. 

o The last Section [7] collects some interesting directions of investigation that merit to 
be developed. For instance the case of Cheeger-type costs in which the product of 
two quantities, scaling in a different way, has to be minimized, as well as the study 
of spectral flows, evolution patterns that follow the gradient flow of a spectral cost 
functional. The case of boundary conditions different from Dirichlet type is briefly 
addressed and referred to the existing literature. 



2 Preliminary tools 
2.1 Capacity 

One of the key tools for treating optimization problems governed by elliptic equations is 
the notion of capacity. 

Definition 2.1. Given a subset E ofW^, we define the capacity of E as 
Cap(S) = inf| [ iVnp + n^dx : uGUe], 

where IAe is the set of all functions u of the Sobolev space H'^{]R.'^) such that u > 1 almost 
everywhere in a neighborhood of E. For every bounded open set D of'K'^ we also define a 
local capacity as 

cap{E,D) = M^j \Vu\'^dx : ueUE^, 

where now Ue is the set of all functions u of the Sobolev space Hq{D) such that u > 1 
almost everywhere in a neighborhood of E. 

Since we are mainly interested in sets with capacity zero the two notions above play the 
same role and in most of the situations they will be equivalent. 

If a property P(x) holds for all x G E except for the elements of a set Z C E with 
Cap(Z) = 0, we say that P(x) holds quasi- everywhere (shortly q.e.) on E, whereas the 
expression almost everywhere (shortly a.e.) refers, as usual, to the Lebesgue measure. 

A subset A of M*^ is said to be quasi- open if for every e > there exists an open subset 

of MJ^, such that Cap{AgAA) < e, where A denotes the symmetric difference of sets. 
Actually, in the definition above we can additionally require that A C A^. Similarly, we 
define quasi-closed sets. The class of all quasi-open subsets of a given set D will be denoted 
hy A{D). 

A function / : D — )• M is said to be quasi- continuous (resp. quasi-lower semicontinuous) 
if for every e > there exists a continuous (resp. lower semicontinuous) function : L) — )■ M 
such that Cap({/ 7^ /e}) < e. It is well known (see for instance [88]) that every function 
u G H^{D) has a quasi-continuous representative u, which is uniquely defined up to a set 
of capacity zero, and given by 

u{x) = lini / / u{y) dy , 

£^0 \B{x,e)\ JB{x,e) 

where B(x,e) denotes the ball of radius e centered at x. By an abuse of notation we 
always identify the function u with its quasi-continuous representative u, so that a pointwise 
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condition can be imposed on u{x) for quasi-every x & D. In this way, we have for every 
E C D 

cap{E, D) = min | j |Vnp dx : u £ Hq{D), u>1 q.e. on 

By the identification above, the quasi-open sets can be characterized as the sets of positivity 
of functions in H^(M.'^); more precisely, this means that A C M"^ is quasi-open if and only if 
there exists a function u € H^(]R.'^) such that A = {u > 0}. 

For every quasi-open set C D we denote by Hq(Q) the space of all functions u G 
Hq{D) such that tt = q.e. on D \Q, with the Hilbert space structure inherited from 
H^iD), i.e. 

Note that Hq{Q) is a closed subspace of Hq{D) (see for instance [76l[88]). Moreover, if Q 
is an open set, the previous definition of Hq{Q,) is equivalent to the usual one (see [2]). 

Most of the well-known properties of Sobolev functions on open sets extend to quasi- 
open sets; for instance, if i}i,i}2 are two quasi-open sets disjoint in capacity, i.e. with 
Cap(Qinfi2) = 0, then H^{niUn2) = H^{ni)nH^{n2) in the sense that u e H^{niUn2) 
if and only if u\q^ £ H^i^i) and Ui^^ G H^i^2)- 

Remark 2.2. All the definitions above have a natural extension to the Sobolev spaces 
VFq'^(O) with 1 < p < +00. We refer to [65] for a review of the main definitions and 
properties of the p-capacity; in particular, when p > d, the p-capacity of a point is strictly 
positive and every '^-function has a continuous representative and therefore, a property 
which holds p-quasi-everywhere, with p > d, holds in fact everywhere, which makes trivial 
several shape optimization problems. This is why, from the shape optimization point of 
view, the most interesting case is when p G {l,d]. 

2.2 F-convergence for sequences of functionals 

In this section we recall briefly the definition and the main properties of the F-convergence, 
which revealed to be one of the key tools for handling variational problems. We do not 
want here to enter into the details of that theory, but only to use it in order to characterize 
the relaxed form of Dirichlet problems; we refer for all details to the books \20\ I52j. 

Definition 2.3. Given a sequence (F„) of functionals from a separable metric space X into 
M we say that (Fn) T -converges to a functional F if for every x £ X 

1. Wxn — T- X F{x) < liminf F„(2;„); 

n— >oo 

2. 3xn —5- X F{x) > limsupF„(xn). 

Theorem 2.4. Properties of T -convergence. 

1. Every T-limit is lower semicontinuous on X; 

2. if {Fn) is equi-coercive on X , that is for every t G M the set IJni-^" — ^'^ relatively 
compact in X, and (Fn) T-converges to F, then F is coercive too, and so it admits a 
minimum on X; 

3. if Xn G argminF„ and Xn ^ x in X, then x € argminF; 

^. from every sequence {Fn) of functionals on X it is possible to extract a subsequence 
V -converging to a functional F on X; 
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5. we have riim„(G + Fn) 
on X. 



G + T lim„ Fn for every functional G which is continuous 



Here we used the notation argmin F to indicate the set of all minimum points of F on X. 
The F-convergence provides a convergence structure on the class T of all lower semicontin- 
uous functionals on X; however, in several situations the functionals under considerations 
fulfill an equicoercivity condition that allows us to introduce a restricted subclass: 



F^{X) = {F:X 



F lower semicontinuous, -F > 



where ^' is a given functional on X. The following result concerning metrizability of the 
F-convergence will be used (see |52j, Theorem 10.22). 

Theorem 2.5. 7/^ : X — )■ M is coercive, i.e. its suhlevels < t} are relatively compact, 
then the class J-\i,{X), endowed with the T -convergence structure, is a compact metric space, 
in the sense that there exists a compact distance 5r on J-iii{X) such that 



5r{Fn,F)^0 



F lim Fn 



F. 



2.3 7-convergence for sequences of quasi-open sets 

Given a bounded domain D C and a right-hand side / G L'^{D) we are interested to 
study how the solution uq, of the elliptic problem 

-Aii = /inJ^, ueHl{VL) 

depends on the domain Vt G A{D). The solutions are extended by zero on D\^}. Since the 
equation above can be also written as the minimum problem 



mm 



D 



iVnl 



f{x)u 



dx : u £ H. 



by the definition of F-convergence introduced in Section 2.2 the study of the limits of (uq^) 
for a given sequence (Qn) of domains can be reduced to the study of the F-limit of the 
sequence of functionals 



rl, 



D 



IVnl 



f{x)u 



dx + In (u) 



where 



ifueH^in) 
+00 otherwise. 



All the functionals above are defined on the separable metric space Hq{D) endowed with 
the L'^{D) convergence. Since the term f{x)udx is continuous, by Theorem 
enough to study the F-limit of the functionals 



2.4 



it is 



Fniu) 



\Vu\^dx + In„{u). 



(2.1) 



D 



Of course, when (0^) converges to in a rather strong way, for instance when r2„ = 
{Id + enV){n) with e„ and F G Cl{R'^;W^), the F-hmit of the sequence Fn above is 
the functional 



Fin) 



D 



\Vu\^ dx + In{u) 



However, in shape optimization problems, we are interested to work with weaker conver- 
gences, in order to ensure, via the direct methods of the calculus of variations, the existence 
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of an optimal domain. Therefore we look for a kind of equilibrium between a convergence 
strong enough to provide continuity, or at least lower semicontinuity, of the cost functionals, 
and on the other hand weak enough to provide the compactness of minimizing sequences. 
When this equilibrium is achieved one obtains the existence of an optimal solution. 
A very natural definition is the following one. 

Definition 2.6. We say that a sequence of quasi-open sets (r2„) in A{D) ^-converges to a 
quasi-open set $1 G -^{D) if the corresponding functionals Jj^ \Vu\'^ dx + Iq,„{u) T-converge 
to the functional |Vup(ix + In{u), in the metric space Hq{D) endowed with the L'^{D) 
convergence. 



By the properties of the F-convergence seen in Section |2.2| we have Q.n — )• O in the 
7-convergence if and only if for every right-hand side / G L'^{D) the solutions uq^^j of 

-Au = fm^}n, U£HQ{il.n) 

(extended by zero on D \ fin) converge in L'^{D) to the solution uqj of 

-Au = fmn, u£H^{n). 

Note that the sequence {uqj) is bounded in Hq{D), hence compact in L^{D). 
Remark 2.7. The following facts for the 7-convergence hold. 



By its definition, and by the equicoercivity of the functionals (2.1 ), the 7-convergence 
is metrizable on A{D). 

It can be proven (see [26j) that the L'^{D) convergence of uq^j to uqj for every / 
holds if and only if it holds for / = 1. Therefore a distance on A{D) equivalent to 
the 7-convergence is 

d^{ni,Cl2) = \\uni,i — uq2,i\\l2(d)- 

For every G A{D) we may define the resolvent operator TZ^ which associates to 
every / G L^{D) the solution uqj; in this way the 7-convergence of (fi^) to Q coincides 
with the pointwise L'^{D) convergence of resolvent operators TZn^ to 

It can actually be proven (see |26j ) that if il„ — >• in the 7-convergence, the con- 
vergence of ur2„,/ to uq^j is in fact strong in Hq{D) and the resolvent operators TZn^ 
converge to TZq in the operator norm C[L'^{D)Y In particular, the spectrum of TZn^ 
converges (componentwise) to the spectrum of TZq, hence the spectrum of —A on 
HQ{Qn) converges (componentwise) to the spectrum of —A on ffg(fi). 



2.4 Capacitary measures 

Unfortunately, the 7-convergence introduced in Section [2. 3| is not compact; in other words 
it is possible to construct a sequence (On) of domains (even smooth) such that the corre- 
sponding solutions do not converge to a function of the form Mn,i- 

The first example of a sequence (On) of this form was provided by Cioranescu and Murat 
in [U]: they constructed the sequence (On) by removing from D a periodic array of balls 
of equal radius Vn 0. If the radius is suitably chosen they proved (see also [26] for 
an interpretation in terms of 7-convergence) that the weak Hq(D) limit of the sequence of 
solutions uq^^i satisfies the PDE 

-An + cti = 1 in Z?, u e H^{D) 
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Figure 1: The sets f7„ in the Cioranescu and Murat example. 



for a suitable constant c > 0, and thus the sequence of domains (iln) cannot 7-converge to 
any domain Q. The sets il„ are represented in Figure [T] 

If we indicate by e„ the distance between two adjacent centers, then the critical size of 
the radius r„ is 

rn = ce^/^"'"^) if d > 2, r„ = exp(-ce-^) if d = 2 

with c > 0. In this way the capacity of the array of removed balls is of order 0(1), while 
the corresponding total volume vanishes as £n 0. 

It is interesting to notice that, for any sequence (f^n) of domains in D and for any 
/ G L'^{D), the solutions uq^j are bounded in Hq{D), so up to extracting a subsequence, 
they weakly converge to some function w E Hq{D). In order to characterize this function 
w as the solution of some limit PDE we have to identify the compactification of the metric 
space A{D) of all quasi-open subsets of D, endowed with the 7-convergence. 

This identification has been obtained by Dal Maso and Mosco in [53] , where it is proven 
that the compactification of the metric space A{D) endowed with the 7-convergence is the 
set A4o{D) of all nonnegative regular Borel measures /i on D, possibly +00 valued, such 
that 

^(B) = for every Borel set B C D with Cap{B) = 0. 

We stress the fact that the measures /i G M.o{D) do not need to be finite, and may take 
the value +00 even on large parts of D. 

Example 2.8. 1. The Dirac measure 5xq does not belong to M.q{D) because it does not 
vanish on the set {xq} which has capacity zero when d > 2. Similarly, the Hausdorff 
measures 'H'^[S, with S manifold of dimension k, do not belong to A4o{D) when 
d>k + 2. 

2. lid — 2<a<d the a-dimensional Hausdorff measure Ti" belongs to Mq{D), and 
consequently every absolutely continuous with respect to Ti" as well. In fact every 
Borel set with capacity zero has an Hausdorff dimension which is less than or equal 
to d- 2. 

3. For every S C D, the measure 005 defined by 

oosiB) = {' ifCap(i?n5) = 0, 
[ +00 otherwise. 
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1/2 



belongs to the class Aio{D). 

Given G AiQ{D) we can consider the relaxed form of the Dirichlet problem by in- 
troducing the space Xf^{D) as the vector space of all functions u £ Hq(D) such that 
JjyU^ dfi < oo. Note that, since /u vanishes on all sets with capacity zero and since Sobolev 
functions are defined up to sets of capacity zero, the definition of X^{D) is well posed. We 
may think to X^{D) as to the Hilbert space (see |38]) Hl{D) n L'^{D, /x), endowed with the 
norm 

ll^^llx^(D) = ( l^^l^ '^^^ Id^^ '^^^ 
which comes from the scalar product 

{u,v)x (D) = / yySIvdx+ I uv dfi. 
Jd Jd 

Given / G L'^(D) we can now consider the relaxed Dirichlet problem 

-An + fiu = f in D, u£ X^{D) 
whose precise meaning has to be given in the weak form 

ueX^{D), / VuVvdx+ / uvdfi = / fvdx \/v ^ X^{D). 
Jd Jd Jd 

The usual Lax-Milgram method in the Hilbert space X^{D) provides, for every // G M.q[D) 
and every / G L'^{D), a unique solution u^j. In this way we can construct the resolvent 
operator TZ^ which associates to every / G L'^{D) the solution u^j. 

Example 2.9. 1. Take n = a{x) dx where a G L^{D) with p > d/2 (any p > 1 if d = 
2). Then, by the Sobolev embedding theorem and Holder inequality, we have that 
Xfi{D) = Hq{D) with equivalent norms. 

2. If is a quasi-open subset of D and fi = oo£)\q then the space X^{D) coincides 
with the Sobolev space Hq{Q,) and the solution u^j defined above coincides with the 
solution uqj of the Dirichlet problem in i}. Thus we can identify the domain to 
the capacitary measure fi = oo£,\q. 

3. If ^ G A^o(-D) and / > 0, then by maximum principle the solution u = TZ^{f) is 
nonnegative, and then also / -|- Au = fj,u is nonnegative. On the other hand, if / > 
we can formally write 

/ + An 

M = 

u 

which gives once u is known; of course it turns out that /x = -|-oo whenever u = 0. 
This argument can be made rigorous (see [46\ I26j ) and therefore, working with the 
class AioiD) is in this case equivalent to work with the class of functions {u G 
Hq{D), u > 0, Au + f > 0}, which is a closed convex subset of the Sobolev space 

The 7-convergence can be extended to the space M.o{D): we have fin ^ fJ- the 7- 
convergence if the solutions 7^^„(/) converge to Tl^{f) weakly in Hq{D) for every / G 
L^{D). Again, it can be proven that this is equivalent to require the convergence only for 

Proposition 2.10. The properties below for the -convergence on the space A4.q{D) hold. 
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1. The space j\4o{D) endowed with the 'y -convergence is a compact metric space; a dis- 
tance equivalent to the ^-convergence is 

d-y{fJ^l,f^2) = -■7^^2(1)||L2(£,). 

2. The class A{D) is included in M.o{D) via the identification Q i— )• oo^i^q and A{D) is 
dense in Mq(D) for the ^-convergence. Actually also the class of all smooth domains 
Q is dense in A4q(D). 

3. The measures of the form a{x)dx with a G L^{D) belong to M.q{D) and are dense 
in Mq[D) for the ^-convergence. Actually also the class of measures a{x) dx with a 
smooth is dense in Mq{D). 

4. If fj-n ^ fJ- for the ^-convergence, then the spectrum of the compact resolvent oper- 
ator TZ^^ converges to the spectrum of IZ^; in other words, the eigenvalues of the 
Schrddinger-like operator —A + fi^ defined on Hq{D) converge to the corresponding 
eigenvalues of the operator — A + ^u. 

Remark 2.11. Even if in this paper we limit ourselves to consider only the case of the Laplace 
operator, the same construction of relaxed domains can be performed for the p-Laplacian 
— Ap in Wq'^{D) for 1 < p < d, and for more general classes of nonlinear operators (see 
[53]). Given / G L'p{D) and a sequence of p-quasi-open subsets of D (i.e. defined in 
a similar way by means of the p-capacity) we denote by Un the solution of the following 
equation on 

-ApUn = / in Qn, Un G WQ^{i}n), 

which has to be understood in the weak sense 

UneWo'^{n„), 1 \VUn\^~^VUnVvdx= f fvdx G Wo'^(f^„). 

A suitable subsequence (still denoted by the same indices) of (un) weakly converges in 
Wq'^{D) to the solution of the equation 

-Apu + fi\u\P-^u = f, n G W^'P{D) n LP^{D), 

H being the Borel measure defined by 

MA)=| ifCap,(^n{» = 0)) = 

\ +00 if Capp(y4 r\{w = 0}) > 0. 

Here w is the weak limit in Wq'^{D) of the solutions above with / = 1 and u = 1 + ApW. 
2.5 Geometrical restrictions 

As we have seen, the 7-convergence makes the space Mq^D) a compact metric space, while 
the subspace of shapes A{D) is not compact (it is even dense in M.o{D)). However, if we 
consider smaller classes than the whole A{D), introducing geometrical constraints on the 
class of admissible shapes, we can provide some extra compactness properties. We refer to 
|26| for all details concerning the restricted classes below. 

• A very strong geometrical constraint is convexity: the class 

A convex of 3-11 convex Open subsets of D 
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is compact with respect to many kinds of convergences, among which the 7-convergence. 



• A weaker geometrical constraint is given by the so-cahed uniform exterior cone condi- 
tion: an open set satisfies the uniform exterior cone condition if for every point xq G dfl 
there exists a closed cone, with uniform height and opening, and with vertex in xq, contained 
in the complement of 0,. The class 

Auni f cone of all Open subscts of D with the uniform exterior cone condition 

is still compact for the 7-convergence. 

• The class Auni f flat cone IS the class of domains satisfying a uniform flat cone condition 
(see [35]), i.e. as above, but with the weaker requirement that the cone may be flat, that 
is of dimension d — 1. The class Auni / flat cone is larger than the previous ones and is still 
compact for the 7-convergence. 

• The previous geometrical constraints can be further weakened by considering the 
classes: 

Acap density of domains Q satisfying a uniform capacitary density condition, i.e. such 
that there exist c, r > with 

cap{Btix)\n,B2t{x)) ^ y,^(n ^ y ^f,o 

.X > c, yte (0, r), Vx G dn 

cap [Bt{x),B2t{x)) 

where Bs{x) denotes the ball of radius s centered at x. 

Auni f Wiener of domains $7 Satisfying a uniform Wiener condition, i.e. such that for 
every point x G dil 

Jr cai.p[Bt{x),B2t[x)) t 
where g : (0, 1) x (0, 1) x D — t- M+ is a flxed function such that for every R G (0, 1) 

lim g(r, R, x) = +oo locally uniformly on x. 

r— >o 

For the classes above the following inclusions can be established: 

Aconvex C Aunif cone C Aunif flat cone C Acap density C Aunif Wiener 

and on each of the classes the 7-convergence is equivalent to the Hausdorff complementary 
convergence Ti^ (i.e. the Hausdorff convergence of the sets D \ Q„) which is then compact. 

• Another interesting class, which is only of topological type and is not contained in any 
of the previous ones, was found by Sverak in [83]; it is concerned only with the case d = 2 
and consists of all open subsets of D for which the number of "holes", i.e. connected 
components of D \ O, is uniformly bounded. More precisely, the classes 

Om{D) = C D, Q open, ^ connected components of D \ < m} 

are such that on them the 7-convergence is equivalent to the Hausdorff complementary 
convergence Ti^, which is then compact. For higher dimensions d> 3 the same result is no 
more true but it can be recovered again if the Laplace operator is replaced by the p-Laplace 
operator, with p £]d — l,d] (see [S])- For p > d the compactness result is trivial, due 
to the fact that points have a strictly positive p-capacity and that the Sobolev embedding 
theorem gives uniform continuity of Sobolev maps. 
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3 Optimization problems for integral functionals 



3.1 Nonexistence of optimal domains 

We present an explicit example of optimization problem for an integral functional where 
the existence of an optimal domain does not occur. We fix a bounded open subset D in M'^ 
(d > 2), a function / S L'^{D), and we consider the minimization problem 

min|j(17) = j j{x,unix))dx : Q £ A{D)^ (3.1) 

where un is the solution of the Dirichlet problem in 0, 

-Au = fmn, ueH^in), (3.2) 

extended by zero on D\^}. 

For instance, taking j{x, s) = |s — uo(x)p, the problem amounts to find a state function 
u as close as possible, in the L?'{D) norm, to the desired state uq, only acting on the shape 



of the Dirichlet region D\Vl, which can then be seen as a control variable, with (3.2) as 
the corresponding state equation. Considering the state function u as the temperature of 
a conducting medium D with prescribed heat sources / in a stationary configuration, the 



optimization problem (3.1 ) consists in finding an optimal distribution of the Dirichlet region 
D\Q.ia order to achieve a temperature as close as possible to the desired temperature uq . 

We shall see that for a large class of integrands j{x,s) an optimal solution does not 
exist; the reason is that a Dirichlet region D\Q composed of many small pieces is more 
efficient, in terms of the given cost, than a domain composed by a single piece. Therefore 
a minimizing sequence (0^) of domains is such that D \ f]„ splits more and more, then 
leading to the nonexistence of the optimum. 



We argue by contradiction and assume that the optimization problem (3.1) admits a 
solution Q that does not coincide with the entire set D; assume also for simplicity some 
mild regularity on $7, as that O is an open set whose boundary has zero Lebesgue measure 
and whose closure O does not fill all the set D. All these extra assumptions could be 
actually removed (see [371 EHl US]) only requiring that is a quasi-open set. For every 
point xo G D \ Q we perform the so-called topological derivative on the cost functional, 
which consists in adding to O a small ball of radius e centered at xq, then obtaining a 
competing domain = il.L) B^. Since for small £ the ball Be is disjoint from Q we obtain 

J{Qe) = / j{x,un{x)) dx + / j{x,Ue{x)) dx + / j{x,0)dx 
Jn JBe Jd\q.i; 

where we denoted by the solution of 

-Au = /inS„ ueHl{Be). 
The same computation for J[Q) gives 

J(0) = / i{x,un{x)) dx + / j{x,Qi)dx+ / j{x,Qi)dx, 
Jn JBe J D\ne 

so that 

J{ne)-J{^)= [ j{x,Ueix)) -j{x,0)dx. 
J Be 

Assuming that j{x,s) is continuous in x and continuously differentiable in s, and that 
the source / is continuous (these conditions can be easily weakened), we obtain 

ue{x) = /(xo) ^'"'^"''"'' + o{e') Vx G Be 
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so that 



j{x,Ue{x)) = j{x,0) +Ue{x)js{x,0) +o{e'^) 

= Jix, 0) + f{xo)js{x, 0) -\x- Xo\^)/{2d) + 0(^2) 
= j{x,0) + f{xo)js{xo,0){e^ - \x-xo\^)/{2d)+oie^). 



Therefore 



Jm - J{n) = I {e' -\x- xo\') dx + o(e'^+2) 



and, using the optimaUty of 17 we end up with the necessary condition of optimahty: 

/(xo)j.(xo,0) >0 VxoeD\17. (3.3) 

For instance, if j{x,s) = \s — uq{x)\'^ and f{x) > on D, the necessary condition above 
becomes 

no(x) < onD\n. 

Taking the desired state uo{x) > on D, we deduce that the only possibihty for a domain 
to be optimal is that = D. However, also this possibility can be excluded if the 
function uq is small enough; in fact, in this case the empty set gives a better value of the 
cost functional, since for small uo we have 

/ liiopdx = J(0) < J(D) = / \uD-uo\^dx. 
Jd Jd 

3.2 Necessary conditions of optimality 

Following the argument illustrated in the previous section we can obtain the result below 
(see ^3138] for a more detailed proof). 

Theorem 3.1. Let f e L^iD) and assume j : D x M — t- M satisfies the conditions: 
1. j{x,s) is measurable in x and continuous in s; 
2- s)\ < a{x) + for suitable a £ L^{D) and b G M; 

3. j{x,-) is continuously differentiable and \js{x,s)\ < ai{x) + 6i|s| for suitable ai G 
L2(J7) and bi G M. 



Then, if a smooth domain ft solves the shape optimization problem (3.1), then the following 
necessary condition of optimality holds: 

f{x)jsix,0) >0 a.e. onD\n. (3.4) 

Remark 3.2. When the optimal set fl is not smooth but only quasi-open, a necessary 
condition of optimality of the form of (3.4) has to be given in terms of the fine topology. We 
recall that the fine topology on D is the weakest topology on D for which all superharmonic 
functions are continuous (see for instance [57] for a detailed study of properties of the fine 
topology); by d*fl and cl*fl we denote respectively the fine boundary and the fine closure of 
fl in D. Under the same conditions on / and j as in Theorem 3.1 the topological derivative 
argument gives (see [371 EE] ) that an optimal quasi-open set Q, must fulfill the necessary 
condition of optimality: 

f{x)js{x, 0) > a.e. on D \ d*Q. 
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Other kinds of necessary conditions of optimality can be obtained for the optimization 
problem (3.1). Assuming again that a smooth optimal domain ^opt exists, it is convenient 
to introduce the adjoint state equation, that for any admissible 0, reads 

— Av = js{x,un) in v G Hq{Q). (3.5) 

In [38] the following optimality conditions have been obtained: 

lit; < q.e. on Qopt', 

f{x)js{x,0) > a.e. on D\ i^opu (3.6) 

{du/dn){dv/dn) = n^-^-a.e. on D n dnopt- 

3.3 Relaxed Dirichlet problems 



We have seen that in general the minimization problem (3.1) may have no solution and its 
relaxed formulation on AiQ{D) 



mm ■ 



.[j{li)= jj{x,u^{x))dx : /iGXo(^)}, (3.7) 

where is the solution of the relaxed Dirichlet problem 

-Au + nu = finD, ueXfj,{D), 

has to be considered, in order to obtain a relaxed solution, which is a measure of the 
class Aio{D). The previous analysis leading to necessary conditions of optimality can be 
performed for relaxed solutions too, considering the adjoint relaxed state equation 

-Av + nv = js{x,u^)mD, v£X^{D). (3.8) 

Let now fiopt an optimal relaxed solution in the space AiQ{D) (which always exists thanks 
to the compactness of the 7-convergence in AioiD)). The optimality conditions (3.6) can 
be written for a general fiopt as 

uv < q.e. in D; 

uv = fi-a.e. in D; (3.9) 

f{x)js{x,0) > a.e. on mt*{u^^^, = 0}, 

where int* denotes the interior with respect to the fine topology. The last condition in 
(3.6), involving normal derivatives, has a more involved translation in terms of capacitary 



measures, and we refer to 

Consider now again the problem of finding the temperature closest to the desired one 
uq, acting on the Dirichlet region D \ 0. We have seen that when uq is small enough and 
the heat source / is positive, no optimal domain exists and the problem has to be written 
in the relaxed formulation 



mm 



\ / \Ufi — uo{x)\ dx : 



in order to obtain a solution in Mq{D). Thanks to the fact that / > and to the identifi- 
cation seen in Example |2.9| (iii), the relaxed formulation above can be rewritten as 



minjy \u - uo{x)\'^ dx : ue H^{D), u>0, An + / > o| 
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which is now a strictly convex problem on the closed convex subset {u G Hq{D), u > 
0, An + / > 0} of the Sobolev space Hq{D). Therefore there is a unique solution u from 
which we may deduce the optimal measure fi through the formula 

/ + An 
= . 

u 

In particular, if I? is a ball and uq, f are radially symmetric, the optimal solution u and 
so the optimal measure fi are radially symmetric too and an explicit computation can be 
made in polar coordinates. Let us consider the particular case when d = 2, D is the unit 
disk, / = 1, and uq = c is a constant. After some computations we find that for small 
values of the constant c the optimal measure fi has the form 

c 

where the values of Rc and of Kc can be computed explicitely by one-dimensional optimiza- 
tion problems. Note that u = c in Br^. A plot of the corresponding optimal temperature 
u{r) is given below. 



0.04- 




Figure 2: Plot of the optimal temperature n(r) for c small. 



4 Spectral optimization problems 

4.1 A general existence result 



We have seen that in general for a shape optimization problem of the form (1.1), either 



we consider on the admissible class A some additional geometrical constraints of the kind 



listed in Section 2.5, or in general the problem does not admit any optimal domain solution 
and the problem has to be considered in the relaxed class of capacitary measures Mo{D). 
However, for some particular type of cost functionals, the existence of optimal domains can 
be obtained (in this case we say that a classical solution exists) without the addition of any 
geometrical constraint but under the only volume constraint \Vt\ < m. In this section we 
show some interesting classes of cost functionals for which the existence of classical solutions 
holds. 

The problems we consider have the general form 

min{F(J]) : O G ^(D), < m} (4.1) 
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where the admissible class A{D) is as before the class of quasi-open subsets of D. A general 
existence result for optimal domains has been obtained in [39] under the following conditions 
on the cost functional F. 

Theorem 4.1. Let F be a 'j-lower semicontinuous functional which is decreasing with 



respect to the set inclusion. Then the optimization problem (4.1 ) admits at least one classical 
solution, i.e. a domain Q G A{D). 

Remark 4.2. Since the 7 is a rather strong convergence (for instance it implies the con- 



vergence of the spectrum, as we have seen in Proposition 2.10) the assumption that F is 
7-lower semicontinuous is very mild and is always fulfilled in spectral optimization prob- 
lems. On the contrary, the monotonicity assumption on F is quite severe; for instance it 
does not hold in the case 

F{Vl)= / \uQ, - uo{x)\^ dx 

JD 

considered in Section [3l 

Remark 4.3. By the monotonicity assumption on the cost functional F in Theorem |4.1[ 



it is clear that the inequality constraint < m in (4.1) can be replaced by the equality 
constraint \^}\ = m. 

Example 4.4. (Optimal domains for integral functionals). Consider a given / G L^{D), 
with / > 0, and let g' : D x M — )•] — 00, +00] be a Borel function such that g{x, •) is lower 
semicontinuous and decreasing on M for a.e. x ^ D, and g{x, s) > —a{x) — bs^ for a suitable 
a G L^[D) and 6 G M. For every 17 G A{D) let = Tln{f) be the unique solution of 
—Am = / on Hq{Q,), and let 

-^(f^) = / 9{x,uq{x)) dx. 

JD 

Then F is lower semicontinuous with respect to the 7-convergence and, since we assumed 
/ > 0, the maximum principle and the monotonicity properties of g imply that F is de- 



creasing with respect to set inclusion. Therefore, by Theorem 4.1 (and Remark 4.3) the 
minimum problem 

min { / g{x,u^{x))dx : Q.^A{D), = m| 
JD 

admits at least a solution. 

Example 4.5. (Optimal domains for spectral problems). For every $7 G A{D) let Afc(il) be 
the k*^ eigenvalue of the Dirichlet Laplacian on Hq(Q), with the convention \k{^) = +00 
if Cap(^) = 0. It is well known that all the mappings Vt 1— )• Afc(il) are decreasing with 
respect to set inclusion; moreover they are continuous with respect to the 7-convergence 



(see Proposition 2.10 (iv) and Remark 2.7 (iv)), so that Theorem 4.1 (and Remark 4.3) 
applies and for every /c G N and < m < \D\ we obtain that the minimization problem 

min{Afe(0) : VL ^ A{D), = m] 

admits a (classical) solution. More generally, arguing in the same way, the minimum 

min{^>(A(J7)) : 9.(^A{D), = m] 

is achieved too, where denotes the spectrum of —A on HQ{yi), that is the sequence 
(Afc(r2)), and the function $ : [0, +00]^ — >• [0, +00] is lower semicontinuous and increasing, 
in the sense that 

A^ ^ Afc VA; G N ^ «>(A) < lim inf $(A'') , 
Xk<i2k V/c G N ^ $(A) < ^{12) . 
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Example 4.6. (Domains with minimal capacity). Since Cap{E) is a increasing set function, 
the mapping Q i— )• F{Q) = Cap(D \ 0,) is decreasing with respect to the set inclusion. It 
is not difficult to verify that the mapping F is also 7-continuous, so that the existence 



Theorem 4.1 applies, and the minimum 

mm{F{n) : n £ A{D), \n\ < m} 

is achieved. If denotes the class of all quasi-closed subsets of D, passing to complements 
and taking m = \D\ — k in the previous problem, we have that the minimum 

min{Cap(£;) : E £ \E\ = k} (4.2) 

is achieved. If we denote by Eq a solution to ( |4.2[ ) we show that 

Cap(^o) = min { Cap(-E) : E C D, |^| = A:}. (4.3) 

In fact, for every subset E of D there exists a quasi-closed set E' such that E C E' and 
Cap(S) = Cap{E') (see for instance Section 2 of |60], or Proposition 1.9 of [5l])- If l-E'l = k, 
then l^'l > A;, so that there exists E" G T with E" C E' and \E"\ = k. By ( |0| )) we then 
have 

Cap(^o) < Cap(^") < Cap(^') = Cap{E), 



which proves (4.3). 



The idea of the proof of Theorem |4.1| is to consider a relaxed solution, that always 
exists, which is a measure /i of the class M.o{D). If we denote by F the relaxed functional, 
defined on A4oiD) we then have 

F(/i) = inf {F(J]) : n £ A{D), \n\<m}. 

We associate to the measure the solution of the relaxed Dirichlet problem 

-Au + fiu = 1 in D, u £ X^(-D) 

and the domain Q = {u^ > 0} which is quasi-open, since £ Hq{D). The new measure 
v = oo£,\fi fulfills the inequality v < ijl and the monotonicity assumption on the cost 
functional gives 

F{n) = F{v)<F{^i), 

then showing the optimality of the domain Q.. 
4.2 The weak 7-convergence 



To give the proof of Theorem 4.1 in a more rigorous way we introduce a new convergence, 
much weaker than 7, that makes the class A{D) compact. We call weak 7 this new con- 
vergence and we denote it by w^. 

Definition 4.7. We say that a sequence (il„) of domains in A{D) weakly ^-converges to a 
domain Q. £ A{D) if the solutions wq^ = TIq^{1) converges weakly in Hq{D) to a function 
w £ Hl{D) (that we may take quasi- continuous) such that Vl = {w > 0}. 



Remark 4.8. We stress the fact that, in general, the function w in Definition 4.7 does not 
coincide with the solution wq, = 7^n(l); this happens only if Vtn 7-converges to il, which 
in general does not occur, because 7-convergence is not compact on A{D). We only have 
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that Qn 7 converges to some G Aio{D), so that the function w in Definition 4.7 coincides 
with the solution tu^ = TZ^{1) of the relaxed Dirichlet problem 

-Au + fiu = 1 in D, ueXf,{D). 

Also, we notice that, by its definition, the i(;7-convergence is compact, since the sequence 
wn„ = T^n„{^) is bounded in Hq{D) so it always has a subsequence (rinj.) weakly converging 
to some function w G Hq(D), and the set of positivity 17 = {w > 0} (which is quasi-open 
since w E Hq{D)) is then the t(;7-limit of (f^Wfe)- 

Finally, by Definitions 2.6 and 4.7 we obtain that the ti;7-convergence is weaker than the 
7-convergence. 

Since the u)7-convergence is rather weak, the class of wj-lowei semicontinuous func- 
tional is much smaller than the class of 7-lower semicontinuous functionals. However, the 
proposition below shows that some relevant examples are still valid. 

Proposition 4.9. Let f £ L^{D) he a nonnegative function. Then the mapping Q 1— )• 
f dx is w^-lower semicontinuous on A{D). 

Proof. Let (O^) be a sequence in A{D) that i(;7-converges to some O G A.{D)] this means 
that wq^ — )• in L^(M^) and that Q = {w > 0}. Passing to a subsequence we may assume 
that Wfi^^ — )• w a.e. on D. Suppose x G is a point where Wfi^{x) — )■ w{x). Then w{x) > 0, 
and for n large enough we have that Wfi^{x) > 0. Hence x G So we have shown that 

^n{x) < liminf In^(x) for a.e. x £ D. 

Fatou's lemma now completes the proof. □ 

In order to show the w^-lowev semicontinuity of other shape functionals, the following 
lemma is needed. 

Lemma 4.10. Let be a sequence of quasi- open sets wy-converging to a quasi-open set 
ri. Then there exists a subsequence (still denoted by the same indices) and a sequence of 
quasi-open sets Gn C D with r2„ C G„ such that Gn ^-converges to $7. 

Proof. Let us denote by wq, the solution 7^n(l) of the Dirichlet problem 

-Au = 1 in 17, u G HI{9)] 

we have w = wq = on D \ 0, and (see for instance [53] for a detailed proof) w < wn on 
ri, and hence w < on the entire set D. 

For each e > we define the quasi-open set fi^ = {wfi > e}. For a subsequence, we may 
suppose that 

wnnUQ" converge to some weakly in Hq (D) 

and, since by a comparison principle we have wq^uq^ ^ wq^, passing to the limit as n — )• 
-|-oo we have that w'^ > wne. Let us show that u;^ G Hq{^}). Indeed, defining = 
1 — ^ min{wQ, e} we get < < 1, u"^ = on 17^, and v'^ = 1 on D \ Q. 
Taking n„ = m[n{v^ , wn^uQe} we get u„ = on $7^ U (D \ (il„ U 17^)), and in particular 
on D\Qn- Moreover Un converges to niin{v^ ,w^} weakly in Hq{D) and hence m.in{v'' ,w^} 
vanishes q.e. on {w = 0}. Since = 1 on D\Q we get that = q.e. on D\Q. Using 
Theorem 5.1 of |54j . from the fact that —Awn^^Qe < 1 in L> we get —Aw^ < 1 and hence 
< wq. Finally wqs <w^< w^, and by a diagonal extraction procedure we get that 

converge to weakly in Hq{D). 

Therefore the quasi-open sets Gn = OnUJl^" 7-converge to 17, which concludes the proof. □ 
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For monotone decreasing shape functionals F the 7-lower semicontinuity and the w'~^- 
lower semicontinuity coincide, as the following result shows. 

Proposition 4.11. Let F : A{D) — )• [—00, +00] he a j-lower semicontinuous shape func- 
tional which is monotone decreasing with respect to the set inclusion. Then F is w^-lower 
semicontinuous. 

Proof. Let us consider a sequence of quasi-open sets in A{D) Ti;7-converging to some 
quasi-open set Q G AiD), and let {J^nfe} be a subsequence such that 

lim = liminf 



By Lemma 4.10 there exists a subsequence (which we still denote by {0^^,}) and D 



such that Gnf, 7-converge to 17. The 7-lower semicontinuity of F gives 

F{n) < liminfF(G„J 

and the monotonicity of F gives 

F(G„J<F(0„J. 

Therefore 

F{n) < liminf J < liminf J = liminf 

fc— >+oo A;— 5>+oo n—>-+oo 

which proves the required zi;7-lower semicontinuity. □ 



Proof of Theorem \4:.l\ The proof of Theorem 4.1 is now straightforward. In fact, if (r2„) 
is a minimizing sequence of quasi-open sets for the optimization problem (4.1), by the 
compactness of the ?i;7-convergence we may extract a subsequence (still denoted by the 
same indices) that tt;7-converges to some quasi-open set G A{D). By Proposition 4.9 we 
have 

iril < liminf |r2„| < m 



and by Proposition 4.11 we have 



F{n) < liminfF(17^ 

n— ^+00 



Therefore J7 is an optimal set for the minimum problem (4.1). □ 
Example 4.12. The minimization 

min{Ai(Jl) : Q e A{D), \n\ < m} 



is one of the first examples of spectral optimization problems. By Theorem 4.1 we know 
that an optimal solution exists and, it was conjectured by lord Rayleigh that, for D large 
enough (to contain at least a ball of measure m) the optimal domain is a ball of measure m. 
The first proof of this fact was obtained by Faber and Krahn by symmetrization techniques 
(see \59\ 1741 [75] ). On the contrary, if D does not contain a ball of measure m, the optimal 
sets Oopt have to touch the boundary dD; it has been shown (see [Ml ES]) that such sets 
^lopt are actually open sets and that their free boundary (i.e. the part of d^lopt included in 
D) is smooth. However, the free boundary D Ci d^lopt does not contain any part of spherical 
surface (see [68l [69] and also ^67j Theorem 3.4.1), in the sense that no part of D n di^opt 
locally coincides with a sphere. 

When the bounding box D is convex we expect that the optimal sets O,opt for Ai are convex 
too; however, this result, even if very natural and strongly expected, is not yet available, 
and the question is still open. 
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Example 4.13. The minimization problem 



min{A2(f^) : n e A{D), \n\ < m} 



also verifies the assumptions of the existence Theorem 4.1 ; as proved in |75tl81j. for D large 
enough (to contain at least two disjoint balls of measure m/2 each) the optimal domain is 
the union of two disjoint balls of measure m/2 each. As before, if D is not large enough 
in the sense above, the optimal sets ^lopt have to touch the boundary dD, but in this case, 
even if it seems reasonable to conjecture that the free boundary of optimal sets is regular, 
the regularity question is still open. Actually, the proof that Clopt are open sets is not yet 
available. 

Minimizing X2{^) in the more restricted class G A{D), \Q\ < m, Q, convexj also admits 
an optimal solution i^opt, as it is easy to prove. When D is large enough, since without 
the convexity assumption the solution is given by two equal disjoint balls, a reasonable 
expectation, also supported in [85] by some numerical computations, is that Qopt is a 
stadium, i.e. the convex hull of two equal disjoint balls tangent each other. This conjecture 
has been disproved in |68t I69| where again it has been shown that, even if O,opt is very close 
to a stadium, the nonflat parts do not locally coincide with a sphere. 

When A: > 3 the optimal shapes for the minimization problem 

m.m{Xk{n) : n e A{D), < m}, 



that exist again thanks to the existence Theorem 4.1, are not known. For /c = 3 and D 
large enough the conjecture, still open, is that: 

the optimal domains for A3 are balls if the dimension d is 2 or 3 and the union of three 
equal disjoint balls when d > 4. 

For A; = 4, d = 2 and D large enough the conjecture, also still open is: 

the optimal domains for A4 in dimension d = 2 are the union of two disjoint balls of 
different radius, whose radii are in the ratio {jo,i/ ji,!)^^"^ ~ 0.79, where jo,i and ji^i are the 
first zeroes of the Bessel functions Jq and Ji respectively. 

For A; > 5 one could expect that optimal sets are also made by a suitable array of balls. 
This is false, as shown by Keller and Wolf for d = 2 and k = 13; in fact in |77j numerical 
computation are provided for d = 2 showing that for k > 5 optimal arrays of disks are not 
optimal (see Figure |3]) . 



Remark 4.14. In Theorem 4.1 the assumption that admissible domains are all contained in 
a given bounded domain D is crucial. A similar existence result for D = M"' is not known. 
For instance, the eigenvalue optimization problem 

i'Pk) mill {Afe(O) : ncR'^, \n\ < m} 

is known to have a solution for A; = 1 (a ball of measure m), for k = 2 (two disjoint balls 
of measure m/2 each) and for A; = 3 (a ball conjectured if d = 2, 3 and three equal disjoint 
balls for d > 4). This last existence result was proved in [33j, where more generally it is 
shown that if the minimizers of problem (Vj) exist and are bounded for j = 1, . . . ,k — l, then 
the minimizers of problem (Vk) exist. For instance, a proof of boundedness of minimizers of 
A3(ri) with Q, C M.'^ would imply the existence of optimal domains for X^IQ,) with Q C M.'^. 

Several other explicit solutions for minimization problems involving eigenvalues are 
known or conjectured; we refer the interested reader to [TUl [^Hl IMl EZl EO] for a wider 
presentation of this subject. 
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No 



Optimal union of discs 



Computed sliapes 



10 



oO 

oOO 



oOO 



OOO 



46.125 



64.293 



82.462 



92.250 



110.42 



127.88 



138.37 



154.62 



oO 



46.125 



64.293 



78.47 



88.96 



107.47 



119.9 



133.52 



143.45 



Figure 3: Optimal domains for and optimal! arrays of disks, 3 < k < 10. 



4.3 Spectral problems with perimeter constraint 

We consider in this section the case when a perimeter constraint is imposed on the admis- 
sible domains. Here we use the De Giorgi definition of perimeter (see for instance [7j): a 
measurable set Q is said to have a finite perimeter if its characteristic function 



Inix) 



1 if X G 
otherwise 



belongs to the space BV(K'^) of functions with bounded variation on M^, i.e. functions in 
L^(M'^) whose distributional gradient is a measure with finite total variation. In this case 
the perimeter of il. is defined by 



Per(^^) 



IVI 



n\- 



Since for a sequence of domains equi-bounded perimeter implies strong compactness of 
the characteristic functions, the situation with perimeter constraint on the admissible class 
could seem a priori better than the one with volume constraint. However the following 
difficulties arise. 

• The relaxation of Dirichlet problems, and hence the appearance of capacitary mea- 
sures as limits of minimizing sequences, may occur even with a perimeter constraint 
on the admissible domains. It is enough to look at the Cioranescu and Murat example 



presented in Section 2.4 the domains 7-converging to the Lebesgue measure can be 
taken with perimeter arbitrarily small. 

The perimeter is not ti;7-lower semicontinuous. Indeed, again by a construction like 
the one of Cioranescu and Murat, we may produce a sequence (f^n) of subsets of the 
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ball B2 centered at the origin and of radius 2 made by removing n small holes from 
the ball Bi of radius 1. If r„ denotes the critical radius in Cioranescu and Murat 
example, it is enough to take the radius s„ of the small holes such that 

In this way, the limit set is = i?2 \ -Bi and 

Per(S2) + Per(Si) = Per(f]) > lim Per(Jl„) = Per(S2). 

• The perimeter constraint gives the compactness of minimizing sequences in the 
convergence of characteristic functions, but the 7-limit can be considerably smaller. 
Indeed, the same example as above, with removed balls radius Sn such that 

r„ < Sn < n'^/'^ 

shows that ri„ — )• B2 in the sense, while r2„ — )• i?2 \ ^1 in the 7-convergence sense. 

The link between tt;7-convergence and L^-convergence is given by the following propo- 
sition. 

Proposition 4.15. Let {An) be a sequence of quasi-open sets which wj-converges to a 
quasi-open set A, and assume that there exist measurable sets 0^ such that An C ^n, o,nd 
that (Cln) converges in to a measurable set Then we have \A\Cl\ = 0. 



Proof. By applying Proposition 4.9 with / = I/j^q we obtain 

\A\Q\ < liminf \An \ ^\= liminf \ = 0, 

n—^oo n— >oo 

which concludes the proof. □ 

In order to consider shape optimization problems with perimeter constraints it is con- 
venient to extend the definition of a monotone decreasing functional F defined on A{D) 
also to measurable sets by setting 

F{A) = mi{F{n) : n C Aa.e., Q G A{D)} . (4.4) 

We notice that, since F is monotone decreasing, we have that its value on A{D) is not 
modified by this extension. 

The shape optimization problems we are interested in are of the form 

mi{F{n) : n C D, Fev{n) < L} , (4.5) 

where L is a given positive real number. 

Theorem 4.16. // F is 'j -lower semicontinuous and decreasing with respect to the set 
inclusion, then there exists a finite perimeter set Clopt which solves the variational problem 
(pl). 



Proof. Let {fin) be a minimizing sequence for problem (4.5). Since Per(i7„) < L we may 
extract a subsequence (still denoted by (f^n)) that converges in to a set Q with Per ($7) < 
L. 

By the construction of the extension of the functional F to measurable sets there are 
quasi-open sets ujn C 0„ a.e. such that 

F{iOn) = F{Qn). 
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By the compactness of tt;7-convergence we may assume that is ?i;7-converging to some 



quasi-open set w, and by Proposition 4.15 we have \ il| = 0. Therefore, we have that 



F{n) < F{uj) < hminf = hminf 



Hence solves the variational problem (4.5). □ 



Of course, if F{0,) = Xi{0,) and D is large enough to contain a ball of perimeter L, the 



balls solve the shape optimization problem (4.5) and are the unique minimizers. Indeed, 



symmetrizing J7 both reduces the cost Ai(r2) as well as the perimeter Per(r2). 

The situation is different for the cost F{Q) = A2(f^), which was first considered in |61j . 
We summarize here the results obtained in [29j, to which we refer for the detailed proofs. 

• In the case d = 2 it is easy to see that the optimal set ^opt has to be convex, since 
convexification both reduces the cost A2 as well as the perimeter. 

• In the case d = 2, if Qopt is a minimizer for A2 with perimeter constraint, then X2{^) 
is simple. 

• Using the facts above, in the case d = 2 it is possible to prove the C°° regularity of 
Qopt- Moreover, if U2 denotes the second eigenfunction, with unitary norm, we 
have the necessary condition of optimality 

Per(Oopt) 

where k{x) is the curvature of dQopt at x. 

• Always in the case d = 2 additional necessary conditions of optimality can be proved, 
as: 

- the boundary of Qopt does not contain any segment; 

- the boundary of Clopt does not contain any arc of circle; 

- the boundary of ^lopt contains exactly two points where the curvature vanishes. 

A numerical plot of i^opt in the two-dimensional case is given in Figure [4| The picture 
shows that minimizers in two dimensions should have two axes of symmetry (one of 
these containing the nodal line), but the proof of this fact is not yet available. 

• If d > 3 no regularity results for the optimal domains are available; actually, at present 
we do not even know if optimal domains are open sets. For a similar problem with 
perimeter penalization the regularity of optimal domains has been proved in |12j . 



If (i > 3 optimal domains are not convex (see [7T]); they should have a cylindrical 
symmetry, even if this fact has not yet been proved. 



5 An example of nonmonotone cost functional 

In this section we consider as a particular cost functional the case of a function which only 
depends on the first two eigenvalues of the Dirichlet Laplacian: 



F(0) = cD(Ai(0),A2(J^)). 



By the results of Section [4. 1| (see Theorem 4.1) we know that the optimization problem 



mm 



|^>(Ai(17),A2(J^)) : neA{D), |l^|<m}, (5.1) 
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^-3-2-1 1 S3 4 

Figure 4: Plot of the optimal set for A2 with perimeter constraint, in the case d = 2. 

without additional geometrical constraints on the admissible domains 17, admits a classical 
solution whenever the function $ is increasing in each of its variables. The following very 
natural question then arises: 

what happens when $ does not satisfy the monotonicity condition above? 

In the rest of this section we suppose the bounding box D is large enough to allow all 
the constructions that will be made. It is convenient to introduce the attainable set 

£ = {{x,y) : x = Ai(il), y = X2{^) for some Q G A{D), \n\ < m}, 

so that the optimization problem ( |5.1[ ) can be rewritten as 

min{$(x,y) : {x,y) G £}. (5.2) 

The set £^ is a subset of which verifies the following properties. 

1. Being the coordinates of points in £ the first and second eigenvalues of the Dirichlet 
Laplacian in some domain J7, we have that x > and y > for every {x,y) G E. 

2. Since A2(f7) > Ai(r2) we have that x <y for every (x,y) G £. 

3. Since the ball B of measure ra makes Ai(r2) minimal, we have x > \i{B) for every 
{x,y) G £. 

4. Since the union A of two disjoint balls of measure m/2 each makes \2{'^) minimal, 
we have y > A2(A) for every (x, y) G £. 

5. Taking the domain Q,/t with t > 1 and using the fact that Xk{il/t) = t^Afc(r2), we 
have that the set £ is conical, i.e. if (x,y) G £ then {tx,ty) G £ for all t>l. 

6. By a result obtained in [11] (proving a conjecture by Payne, Polya and Weinberger 
stated in [78]) the balls minimize the ratio Ai(r2)/A2(f7) among all domains Q, there- 
fore y < xX2{B) / Xi{B) for every {x,y) G £. 
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A numerical output of the set £ in the two-dimensional case has been obtained in [87] and 
is reported below. 




Figure 5: A numerical plot of the set £. 
If the function $ is lower semicontinuous on and satisfies the coercivity condition 

lim ^{x^y) = +00 (5-3) 

\{x,y)\'^+oo 



then the existence of a solution to problem (5.2), and then to problem (5.1), follows straight- 
forward by the Weierstrass theorem as soon as we can prove that the set £ is closed in R^. 

Proposition 5.1. // the set E is convex, then it has to be closed. 
Proof. Let (x, y) £ £ and consider the minimization problem 

min {(x — x)"*" + (y — y)^ : (x,y)G<5}, 

which is equivalent to 

mm {{Xi{n) - x)^ + {X2in) - y)^ : Q £ AiD), \n\<m}. 
Since the function (x — x)"*" + (y — y)~^ is increasing in each of its variables, by the existence 



Theorem 4.1 the minimum problems above ctdmit optimal solutions {x opt "^y opt 

) G £ and 

^opt £ -^{D) respectively. Since (x, y) G <5 it is also clear that the minimum value has to 
be zero, so that 

Xopt ^ X, Vopt ^ V' 

If we are assuming that the set £ is convex, the segment joining {xopt,yopt) to the point of 
£ given by [Xi{B), X2{B)), where B is the ball of measure m, is all contained in £, as well 
as the segment joining {xopt,yopt) to the point of £ given by [Xi{B), X2{B)) , where B is the 
union of two disjoint balls of measure m/2 each. By the conicity property (v) above of the 
set £ the point (x, y) has to belong to £ thus proving that £ is closed. □ 
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Unfortunately, in spite of the numerical evidence of the convexity of £ provided by 
Figure [5| a proof of the convexity of £ is still missing, so the result of Proposition 5.1 above 



is useless to deduce the closedness of £ and then the existence of optimal domains for the 



minimization problem (5.1). However, a result stronger than the one of Proposition 5.1 
holds. 

Definition 5.2. We say that the set £ is: 

• horizontally convex if for every point (xo,yo) £ £■ the horizontal segment joining 
(xotUo) to the straight line y = x is all contained in £; 

• vertically convex if for every point (xo,yo) ^ ^ the vertical segment joining (xo,yo) to 
the straight line y = kx, with k = X2{B) / \i{B) being B any ball of measure m, is all 
contained in £ . 

Of course, if £ is convex it is also horizontally convex and vertically convex, while the 
converse could not in principle be true. It happens that horizontal and vertical convexity 
are still enough to prove that £ is closed. 

Proposition 5.3. // the set £ is both horizontally convex and vertically convex, then it has 
to be closed. 



Proof. It is enough to repeat the proof of Proposition 5.1 and to conclude that {x,y) S £ 



by using horizontal and vertical convexity instead of convexity. □ 

While the convexity of the set £ is still unproved, the horizontal and vertical convexity 
can be showed; we give a sketch of the proof omitting some technical details. A complete 
proof can be found in f27j . 

Proof of horizontal convexity. We have to show that, given a domain 0, G A{D) with < 
m, we can construct domains i}{t) S -^{D) with \i^{t)\ < m such that 

\2{n{t)) = \2{^), \i{n{t)) = t\2{^) + {l-t)Xi{^), VtG[0,l]. (5.4) 

Assume that is a regular set; if we denote by a second eigenfunction for Q and by S the 
corresponding nodal set, i.e. the set {n = 0}, we may continuously shrink S, obtaining then 
for t G [0, 1] sets S{t) C S continuously shrinking from S to 0. The sets 0,{t) = Q, \ S{1 — t) 
then 7-continuously move from 0, to ^} \ S. Since S is a nodal set we have 

X2{n{t)) = X2in) vtG[o,i]. 

On the other hand, since 7-continuously decreases from Q to Q \ S, we have that 
Ai(r2(t)) continuously increases from Ai(r2) to Xi{Q \ S) = X2{^). 



A reparametrization of the map 1 1— )■ Q{t) now gives the required property (5.4). □ 



A more careful analysis is made in [27] where the following result is shown. The appli- 
cation to horizontal convexity is then as above, with = Q and J^i = $7 \ 5. 

Proposition 5.4. Let Qi C be two quasi-open sets. Then there exists a decreasing 
homotopy from to Qi which is "y- continuous, namely there exists a 'j- continuous mapping 
h : [0, 1] A{D) such that 

h{t) C h{s) for s < t, h{0) = Qq, h{l) = Qi. 
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An important tool that we shall use in the proof of horizontal convexity is the so-called 
continuous Steiner symmetrization (see |24l I25j). Roughly speaking it consists in a path 
1 1— 7- fit starting at t = from any domain fi £ -4(D) and symmetrizing it more and more, 
finishing at t = 1 in a ball of the same measure. The crucial fact is that during the evolution 
the first eigenvalue Xi{^lt) decreases, as a consequence of the fact that Qt is more and more 
symmetric. 

Proof of horizontal convexity. We have to show that, given a domain 0, G A{D) with \Q\ < 
m, we can construct domains fij G •^(D) with |fit| < m such that 

Ai(fit) = Ai(fi), Ai(fit) = tA2(S) + (l-t)A2(fi), VtG[0,l] (5.5) 

where, again, B denotes the ball of measure m. If the continuous Steiner symmetrization 
fit was 7-continuous, the proof could be easily achieved because the path 

x{t) = Ai(a), y{y) = A2(fit) 

would join continuously the point (Ai(r2), A2(fi)) to the point (Ai(i?), A2(-B)) with x{t) 
descreasing, and this, together with the conicity property (v) of the set £ would allow us 



to conclude that for every t £ [0, 1] there is a set fit verifying the property (5.5). Actually, 
the continuous Steiner symmetrization path fit provides a mapping t i— )• A2(fit) which is 
only lower semicontinuous on the left and upper semicontinuous on the right, but this is 
still enough to achieve the proof, using again the conicity property (v) of the set £. □ 

Summarizing, we have obtained the following existence result 
Theorem 5.5. If D is large enough, and $ is lower semicontinuous and verifies the coer- 



civity condition (5.3), then the shape optimization problem (5.1) admits a solution. 



In the proof, several properties of the first two eigenvalues of the Laplace operator have 
been used; we do not know if more general problems of the form 



mm 



{«>(A,,(fi),...,Ai,(fi)) : neAiD), |fi|<m} 



still admit an optimal domain solution, without monotonicity assumptions on the function 
for different choices of the indices ii, . . . , ik- 

6 Optimal partition problems 
6.1 Existence of optimal partitions 

In this section we consider the minimization problem for shape cost functionals of the form 

F(fii,...fi^) 

where is a fixed integer, the unknown domains fij G A{D), the cost is a map F : 
A{D)^ — )- [0, +oo] and the domains fij fulfill the condition 

fij n fij = for i / j. 

A family {fii,...,fi7v} of pairwise disjoint subsets of D will be called a partition. The 
minimization problem we consider is then 

min{F(fii,...fi7v) : Qj £ A{D), D = Hi for i ^ j} . (6.1) 

As in the case of a single domain fi we may consider several interesting cases of partition 
cost functionals F: 
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1. integral costs given by 

F{n,,...nN)= [ j{x,nnAfi),---,'^nJfN))dx 
Jd 

where fi £ L^{D) are given and as usual the solutions TZuiifi) of 

—An = fi on Qi n G -ffQ(f]j) 
are intended as extended by zero outside fij 

2. spectral costs given by 

F(17i,...f^jv) = ^(A(ni),...,A(f^^)) 
where we denoted by X{Q) the spectrum of the Dirichlet Laplacian in Q. 



By the same methods used in Section |3.2| we may obtain some necessary conditions 
of optimality for the case of integral costs (i) above. Assuming that an optimal partition 
(ill, ■ ■ ■ ,^n) exists and is made of smooth domains we have 

UiVi < q.e. on Qi] 

fi{x)js,{x,0,...,0)>0 a.e.onD\ni; (6.2) 
{dui/dn){dvi/dn) = n^-^-a.e. on L» n {d^i \ IJ^^j 9%) 

for alH = 1, . . . , A^, where we denoted by Ui and Vi the solutions of the state and adjoint 
state equation 

—Am = /j on fij MS HQ{ili) 

and 

-Av = jsiix,ui, . . . ,un) on^li V e i?o(flj) 

respectively. 

In general, without extra conditions either on the admissible partitions or on the cost 
functional F, we could reproduce counterexamples to the existence of optimal partitions 
similar to the ones seen in Section 3.1, and optimal solutions only exist in a relaxed sense. 
For optimal partition problems the relaxed formulation will be considered later in Section 
|6.2[ here we focus our attention on some monotonicity assumptions on the cost F that will 



imply the existence of classical solutions to problem (6.1). 
Definition 6.1. We say that F : A{D)^ [0, +oo] is: 

1. ^-lower semicontinuous if 



F{ni,...,nN) < iimmfF{n'i,...,n% 

n—>-+oo 

whenever $7" — )• in the j -convergence, for i = 1, . . . , N ; 
2. wj-lower semicontinuous if 

F{ni,...,nN) < liminfF(0^,...,0^ 

n—>-+oo 

whenever $7" — )• ilj in the wy -convergence, for i = 1, . . . , N ; 
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3. monotonically decreasing in the sense of the set inclusion if for all {Ai, . . . , A^), 
{Bi, . . . , Bn) G A{D)^ such that Ai C Bi for alii = 1, . . . , N in the sense of capacity, 
i.e. Cap{Ai \ Bi) = 0, then 

F{Bi,...,Bn) <F{Ai,...,An). 

The following existence result for optimal partition problems with monotonically de- 
creasing shape cost functional has been obtained in [28]. 

Theorem 6.2. Let F : A{D)'^ — t- [0,+oo] be a shape cost functional which is wj-lower 



semicontinuous. Then the optimal partition problem (6.1) admits a classical solution. This 
happens for instance if F is j-lower semicontinuous and monotonically decreasing in the 
sense of the set inclusion. 

For instance Theorem |6.2| applies to the case 

F(Oi, ...,nN) = ^'(Afci(Oi), . . . , XkA^N)) 

where ki > 1 are integers and $ is increasing in each variable, i.e. 

<^(si, . . . , sn) < ^{ti, . . ■ , tiy) whenever Si <ti, i = 1, . . . , N. 
Note that no measure constraints on the sets ili are imposed; however, thanks to the 



w7-lower semicontinuity result of Proposition 4.9, if F : A{D) — )• [0, +oo] is tf7-lower 



semicontinuous, the measure constrained minimization problem 

min {F{^li, ...Qn) ■ & -^(.D), ^li n = for i / j, < nii] 

still admits a classical solution. 

A particularly interesting example is given by the shape cost functional 

= Ai(l^i) + --- + Ai(J]iv). 

In [43| m] Caffarelli and Lin considered the equivalent variational formulation of the prob- 
lem 

min I j \Vu\^dx : u€ H^{D;R^), j uf dx = 1 for sdl i, G(n) = o} 

where 

and showed the regularity of the free boundary of the optimal partition, up to a singular set 
of Hausdorff dimension less than or equal to d — 2. They also formulated a very interesting 
conjecture about the behaviour of the optimal partitions when the number of parts tends 
to +oo: 

in the two-dimensional case d = 2 the optimal partitions {$7°^*} tend, as N ^ +00, to be 
made by a regular exagonal tiling; in particular, for the optimal value Pjsi = "^1^=1 ^li^T^) 
the estimate 

holds, where H is a regular exagon of area equal to 1. 

The numerical computations made in [19] confirm the conjecture, as Figure [g] shows. 
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6.2 Relaxed formulation of optimal partition problems 

For a general partition functional F{^}i, . . . which does not satisfy any monotonicity 

assumption we cannnot expect the existence of a classical solution to the optimization 



problem (6.1). Similarly to what happens in the case of a single domain, a relaxation 



procedure is needed to describe the behaviour of minimizing sequences, and this will be 



done through the use of the capacitary measures described in Sections 2A and 3^ 

In order to characterize the expression of the relaxed problem associated to (6.1) we 
consider, for every sequence {^2" , . . . , r^^} of pairwise disjoint quasi-open subsets of D the 
associated measures ^u" = oo^j^^n (i = 1,. . . ,N) of the class J^o{D) introduced in Section 



2.4 Since Mo{D) is compact with respect to the 7-convergence, up to a subsequence, there 
exist measures /ij G M.o{D) (i = 1, . . . , N) such that fii is the 7-limit of /i". 

The fact that the original $7" were pairwise disjoint will imply some conditions on the 
limit measures /Uj, that are not "independent", and characterize the relaxed optimization 



problem associated to (6.1). For instance, it is not possible to obtain an A^-tuple made by 
all the measures /i, equal to the Lebesgue measure on D. 

For every capacitary measure n £ M.o{D) we denote by 0^ the set "of finiteness" of ^, 
more precisely defined as 

^7^ = {7^^(l) >o}. 

The following result has been proved in |41j . 

Theorem 6.3. An N -tuple [fii, . . . , hn) of capacitary measures is made by the ^-limit 
of pairwise disjoint quasi-open sets {17", . . . , il^} if and only if it satisfies the following 
property: 

Cap(0^, n J7^J = Vi/j. 

For instance, if 



f{vii,...vln)= I j(x,7^Ql(/l), . . . ,7^^^(/A^)) dx, 

Jd 

then the relaxed problem associated to 

min . . . ^n) ■ G AiD), ^ n = for i / j} 

is given by 

min{F(^i, . . ./iAf) : £ Mo{D), Cap(Q^^ n fi^^) = for i 7^ j} 

where 

. . ./iAf) = / i{x,n^^{fi),...,n^^{fN))dx. 
Jd 
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7 Further problems 



7.1 Cheeger-type problems 

We consider here minimization problems that we may call of Cheeger type, for the similarity 
with the well-known Cheeger problem 



mm 



It is possible to prove that for every bounded domain D there exists an optimal Cheeger 
set ^opt and (see |31I35]) this set is unique and convex whenever D is convex. Moreover, 
in this case the boundary dflopt does not contain the points of dD with too large mean 
curvature; more precisely, di^opt coincides with dD if and only if 

mil <^(^) 

where H{x) is the mean curvature of dD at x and X{D) is the minimal value of the Cheeger 
minimization problem above. Note that the two quantities Per(r2) and \^}\ scale in a different 
way and, due to the isoperimetric inequality 

Per(f^)|f]|(i-'^)/°' > Per(S)|B|(i-^)/'^ = du;]^'^ 

(where uj^ denotes the Lebesgue measure of the unit ball in M'^), we have 

Per(O) 
hm ^ = +00. 

|n|->o 

The same analysis occurs for the cost 

|r2|~" Per(i7) whenever a > '^—^ — . 

In this section we consider rescaled minimization problems of the form 

min{M(J7)J(0) : ^ e A(D)] (7.1) 

where the mappings M and J fulfill some rather general assumptions related to the varia- 
tional 7-convergence on the family A{D) of quasi-open subsets of D. These problems have 
been studied in [42j to which we refer for all the details. On the mappings M and J, from 
A{D) into [0, -|-oo] we assume 

1. M and J are nonnegative, and J{D) > 0; 

2. J is 7-I.S.C. and decreasing with respect to the set inclusion; 

3. M is w"f-\aweT semicontinuous; 

4. the Cheeger scaling condition is verified: 

lim M{Q)Ji^) = +00. 
M(n)^o 

With this last assumption we may define the cost M{Vl)J{^) = +00 whenever M(il) = 
and obtain the following existence result. 
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Theorem 7.1. Under the assumptions above the minimum problem (7.1) admits a solution 
Qopt and M{0,opt) > 0. 

Proof. If (Sin) is a minimizing sequence, by the compactness of the t(;7-convergence (see 
Section |4.2| ) we may assume that Jin tend in to some SI G A{D), with iVf(r2) > 
by the assumptions above. By the properties of t/;7-convergence and by Proposition |4.11 



the functional J is u;7-lower semicontinuous, as weU as the product MJ, which allows to 
conclude that is an optimal domain for the minimum problem (7.1). □ 



Example 7.2. The following cases fulfill the assumptions of Theorem 7.1 



Let J{Q,) = <I>(A(r2)) where X{Q) denotes the spectrum of the Dirichlet Laplacian in 
and $ is a lower semicontinuous function increasing in each variable and positively 



p- homogeneous. Taking M(r2) = ir^l"^ we have that the assumptions of Theorem 7.1 
are fulfilled whenever a < 2p/d, so that in this case the minimization problem 

min{|ll|"$(A(ll)) : Q e A{D)} 

admits a solution. 

The case J{il.) as above with M(S7) = |Per(r2)|" does not fall into the framework 
of Theorem |7.1| because the shape functional Per(S7) is not wj-lower semicontinuous 



(see Section 4.3). However, using the arguments of Section 4.3 we obtain that the 
minimization problem 

min{|Per(0)|"$(A(0)) : Q e A{D)} 
admits a solution whenever a < 2p/{d — 1). 
Consider the integral functional 

I{n)= [ a{x)\un{x)\P dx 
Jd 

where uq is the solution of the Dirichlet problem in ^1 

-Au = fmn, u£H^{n), 

extended by zero on D\Q. Assume that a G L"^ is nonnegative, with q > p/{p — 1), 
/ G is nonnegative, and < p < 2d/{d — 2). Then the shape functional 

j{n) = i/i{n) 



fulfills the assumptions of Theorem |7.1| and, taking M{i}) = |S7|° we have that the 
minimization problem 

admits a solution whenever a < (2p + d) /d. A similar analysis holds for the shape 
functional 

Jin) = [ a{x)\un{x)\^P dx. 
Jd 
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• The case of the minimization problem 



mm 



where /(O) is as above, does not fall into the framework of Theorem 7.1 but the 
arguments of Section [4.3| apply, and we can conclude that an optimal domain exists, 
provided that a < {2p + d)/{d — 1). 

As it happens in the Cheeger problem, due to the scaling condition we assumed, the 
optimal domains ^opt for the examples above have to touch the boundary dD. Here below 
we describe some necessary conditions of optimality, obtained in |42] for the minimization 
problem 

where — Auq = 1 in Jl, with Dirichlet boundary conditions. 

Theorem 7.3. Assuming that ^lopt is smooth, and that a < (2 + d)/d, the following nec- 
essary conditions of optimality hold: 

\Vun^J^ = aK on dQ^pt n D K - ^ f uo dx 

Remark 7.4. Assume that dD has isolated conical points, but is smooth otherwise. The 
gradient of the solution un^p^ then vanishes in the conical points and a pointwise decay 
estimate for the gradient holds. Theorem 7.3 then shows, that the optimal domain O,opt 



cannot fill the entire set D. 

In a similar way we can obtain necessary conditions of optimality for the minimum 
problem 

min{|J7|"Afc(17) : n G A{D)^ 
where a < 2/d. If we denote by u a normalized k-th eigenfunction of ^lopt we obtain (see 

iVul"^ = aH on dQopt(^ D; _ Xkj^opt) 

iV-up > aH on dQopt n dD. '''''' \n^t\ ' 

We conclude this section by pointing out some shape optimization problems for which 
the existence of a solution is still unavailable. For a fixed k > 1 consider the minimization 
problem 

mm{XU^) I undx : n G A{D)} 
Jd 

with a > 1 + d/2. By the results of [721 E3] we have 

Ai+'^/^(Q) / undx> xI^'^^\b) [ UBdx 

JD JD 



for any ball so that the scaling condition 



lim M(n)J(n) = +00 



is fulfilled whenever a > 1 + d/2, taking 

J{n) = X'^{n) and M{n)= [ Undx. 

JD 

However, the t(;7-lower semicontinuity condition fails for M(r2), as it can be easily seen by 
a sequence of finely perforated domains like the Cioranescu and Murat ones (see Section 



2.4|, and so the existence Theorem 7.1 cannot be applied. It would be interesting to prove 



(or disprove) that an optimal domain flopt for the problem above exists. 
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7.2 Spectral flows 

In [56] De Giorgi introduced a very general theory to study evolution problems with an 
underlying variational structure. The theory was called minimizing movement theory and 
its framework is so flexible to be applied both to quasistatic evolutions as well as to gradient 
flows, under rather mild assumptions. Here we limit ourselves to recall the part of the 
scheme which is interesting for our purposes, referring to [8J for further details. 

Let {X, d) be a complete metric space, let uq G X he an initial condition, and let 
F : X ^] — oo, +00] be a functional defined on X. For every fixed e > the implicit 
Euler scheme of time step £ and initial condition uq consists in constructing a function 
U£{t) = w{[t/£]), where [•] stands for the integer part function, in the following way: 

{d (v w(ti)^ 1 
-^(^) "I '2 /■ 

Definition 7.5. If a is another topology on X, we say that u : [0,T] ^ X is a generalized 
minimizing movement (or a variational flow) associated to F and a, with initial condition 
uq, and we write u £ GMM{F,(t,uq), if there exist a sequence e„ — ^ such that 

Ue„{t) u{t) in a Vt € [0,T]. 

The simplest situation occurs when a is the topology of the metric d; in this case the 
following existence result for variational fiows holds (see Proposition 2.2.3 of [8]). 

Theorem 7.6. If F is d-lower semicontinuous and d-coercive, in the sense that its sublevels 
{F < k} are d-compact in X , then for every initial condition uq £ X there exists u G 
GMM{F,d,UQ). Moreover, u belongs to the space ^C^((0,T);X) of mappings such that, 
for a suitable m G L^(0,T) 

d{u{s),u{t)) < [ m{r)dr Vs,tG[0,r]. 



The result above applies straightforward, thanks to Proposition 2.10 (i), to the case X 
Mq{D) endowed with the compact distance d^. By Theorem 7.6 for every initial condition 



/Uo eMo{D) there exists a GMM{F, d^,fio) flow n{t) which is of class AC'^ ((0, T);Mo{D)) . 

There is a natural one-to-one map between the class of capacitary measures A4o{D) 
and the closed convex set of L'^{D) 

K = {w€H^{D) : u; > 0, 1 + Aw>0}, 

given by 

^ , , , 1 + Aw 

u I—)- Hall) := Wn, with mverse w ^ a-a, = . 

w 

Moreover, the metric structures on Mq{D) and K are the same, since 

djifJ-1,1^2) = - ■^«M2llL2(D)• 
Therefore, every functional F : Aio{D) — )■ M can be identified with a functional J : X — >• M 

by 

F{p) = J{wfj) or equivalently J{w) = F(/i^). 

The variational flow for F in Mq{D) can be then obtained through the variational flow for 
J in the metric space K endowed with the L'^{D) distance, generated by the implicit Euler 
scheme 

G argmin I J(t(;) -h ^ y - tj;^^ dxj. (7.2) 
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For instance, we may consider a functional J of integral type, of the form 



J{w)= / j{x,w{x)) dx, 
Jd 



ID 

where j : D x M — t- M is a suitable integrand. In particular, we may take j{x,w) = —w 
which leads to the energy of the system for the constant force / = 1, oi j{x,w) = w which 
gives the torsional rigidity. If j{x,s) is convex and lower semicontinuous in the second 
variable the variational flow for capacitary measures is reduced to the gradient flow of a 
convex lower semicontinuous map on L^(D). 

An interesting question is the following: if we start from an initial condition which is 
a quasi-open set, i.e. /xq = oo£,\Qg, in which cases the flow remains in the family of quasi- 
open sets? It can be shown that in general this does not happen, at least at the discrete 
level. This kind of phenomenon was numerically observed in the framework of quasi-static 
debonding membranes |30j . 

We may consider, instead of capacitary measures, flows of shapes by endowing the 
class A(D) of quasi-open subsets of D with a suitable distance. There are no "standard" 
distances on A{D) and several choices can be made. A first possibility is to consider the 
Lebesgue measure of the symmetric difference set 

Since two quasi-open sets may differ for a negligible set (as for instance in a disk and 
a disk minus a segment), this is not a proper metric in A{D), so that one should consider 
equivalence classes in the family of shapes. 

The distance dchar is not compact on A{D); nevertheless, for 7-lower semicontinuous 
functionals which are monotone decreasing for the set inclusion, we have the following 
result. 

Theorem 7.7. Let F : A{D) M be a j-lower semicontinuous functional monotone 
decreasing for the set inclusion, and let S A{D). There exists a GMM map t 1— )• 
associated to F and to the distance dchar: with initial condition Qq. Moreover, the flow 
is increasing for the set inclusion. 

For instance, spectral functionals of the form 

F(l^) = cD(A(17)), 

where X{Q) is the spectrum of the Dirichlet Laplacian in and <l> : M'^ — t- M is increasing 
in each variable and lower semicontinuous, fulfill the assumptions above. 
Several questions about the fiow arise. 



The fiow is not obtained through Theorem 7.6 so the functional F could be 
discontinuous on the curve it would be interesting to analyze this issue and to 

identify some cases where this continuity occurs. 

In M^, assume that Qq is simply connected. Prove or disprove that the GMM associ- 



ated to Ai in the framework of Theorem 7.7 consists only on simply connected open 
sets. 

Assume that 0,q is convex. Prove or disprove that a GMM associated to Ai in the 



framework of Theorem 7.7 consists only on convex sets. 
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• Assume that CIq is convex. Is it true that any minimizer of 



min Ai(0) + 1171 



is convex; 



Is it true that the GMM associated to Ai in the framework of Theorem 7.7 converges, 
up to a suitable rescahng, to a ball? 

Prove or disprove that the metric derivative of Ai, computed at a bounded smooth 
set Q is given by {ui denotes the normalized first eigenfunction) 



dui 
dn 



Other possibilities for defining spectral flows are possible, as for instance the frameworks 
below. 

o Work in the class of sets with prescribed measure {Vt G A.{D) : \^\ = m}. 

o Introduce a penalization on the perimeter, that is work with the functional F{^}) + 
FeT{n). 

o Work with the Hausdorff complementary distance (i//c in the family of open subsets 
ofD: 

dH<=i^i, ^2) = max D \ ili) — d{x, D \ Cl2)\- 
xeD 

o Work in the class of convex open subsets Q metrized by one of the distances: 

- The Hausdorff distance; 

- The distance d^har of the characteristic functions; 

- The distance of the oriented distance functions. 



7.3 Other kinds of boundary conditions 

All the results presented in the previous sections are concerned with the framework of 
Dirichlet boundary conditions. The main reason why this framework allows to treat the 
various questions we have seen, is the monotonicity relation 

f>0,nicn2 7^nl(/) <^f72(/) and A(1^2) < A(l^i). 

When working with different boundary conditions, the monotonicity properties above are no 
more true, which creates several difficulties if one wants to obtain similar kinds of results. 
In addition, while a function u G ^o(^^) is easily extended to Hq{D) by putting u = 
outside fi, similar extension operators do not exist if the Dirichlet boundary condition is 
replaced by Neumann or Robin type conditions. Finally, in order to work with solutions of 
elliptic PDEs with Neumann or Robin boundary conditions, some mild regularity on the 
domains O has to be required, which could be lost when passing to the limit on sequences 

Nevertheless, some shape optimization results are known in these more difficult frame- 
works; we mention few of them together with some open problems, referring for instance 
to [67j for more details. 

• Neumann boundary conditions. If we denote by fJ.{^) the spectrum of the Neumann- 
Laplacian on the domain $7 (notice that = for all f]), the minimization problem 
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for under the volume constraint = m has a trivial solution. Indeed, we have 

= for any domain Q which has at least k connected components. 

If we impose that the domains Q must be convex and with a given diameter, then the 
infimum is not zero, but it is not achieved, as shown in |79j. 

A more interesting problem is the maximization of fiki^) under a volume constraint. 
It is known (see \84:\ I86j) that for /i2(f^) the maximum is attained for a ball, which are 
the only maximizers. An analogous result for A; > 3 (i.e. existence of optimal domains 
and possibly their characterization) is not known. Similarly, the existence of a maximizing 
domain for /i2(^^), with the constraint = m, is not known if we impose the constraint 
Q C D with D narrow enough to not contain a ball of measure m. 

Some more issues on spectral optimization problems for the Neumann-Laplacian, to- 
gether with related numerical computations, can be found in [15]. One could also consider 
mixed boundary conditions: Dirichlet on one part and Neumann on the remaining part. We 
mention the paper [49j where this problem has been studied, and [9j where the behaviour 
of the spectrum of the Neumann Laplacian under boundary perturbations is analyzed. Fi- 
nally, some results are available for minimization problems of functions of reciprocals of fik 
under volume constraint (see [67]) as: 

1 1 V ^ V ^ 



• Robin boundary conditions. If is a sufficiently regular open subset of M'^ the Robin 
eigenvalue problem for the Laplace operator is given by 

du 

—An = rju in $7, — — h f3u = on Oil, 
on 

where /3 > is a given real number. It is known (see for instance |55| [STj \T7\ I18| ) that, 
denoting by ??i(r2) the smallest Robin eigenvalue, the minimum of ??i(r2) among the class 
of bounded Lipschitz sets with prescribed volume is achieved on balls, which are also the 
unique minimizers. 

However, the problem of setting spectral optimization problems with Robin boundary 
conditions for general nonsmooth domains, faces some difficulties dues to the presence of 
boundary terms. In [32] a new approach has been proposed, which, using the features of 
SBV functions, introduced by De Giorgi for treating free discontinuity problems (see [7]), 
allows to define tii{Q) also for nonsmooth domains. Denoting by SBV^^'^(R'^) the space 

SBV^^^iR'^) = {u : R'^ ^ [0, +oo[ measurable and n^ G SBViR'^)} 

we may define i]i(Q) for a general open set through the Rayleigh type formula 

77i(J7) = minj /" \Vu\^dx + P f [{u+f + {u-f]d'H'^'^: [ dx = l} 

where the minimum above is taken over the class of functions u S SBV^^'^{M.'^) with | spt ii\ 
r?| = 'H'^~^{Ju \ do,) = 0, and Vu is the approximate gradient of u, Ju is the set of 
discontinuity points of u, are the traces of u on Ju from the two sides. In this way in 
|32j it is proved that for every open set O the isoperimetric inequality 

rim > VliB) 

holds, where i? is a ball with \B\ = Moreover, equality holds if and only if is a ball 
up to a negligible set. 

However, all the questions of general spectral optimization problems still remain open, 

as: 
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- existence results under the constraint Q C D; 

- treating more general spectral costs as $(r/(ri)); 

- regularity of optimal domains. 
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